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Non-Hermitean operators may appear during the calcula- 
tion of a partition function in various models of statistical 
mechanics. The tail eigen-states, having anomalously small 
real part of energy Re{e), became naturally important in this 
case. We consider the distribution of such states and the form 
of eigenfunctions for the particle propagating in an imagi- 
nary random potential (the model motivated by the statis- 
tics of polymer chains). Unlike it is in the Hermitean quan- 
tum mechanics, our tail states are sufficiently extended. Such 
state appear if the values of random potential turns out to be 
anomalously close inside the large area. Results of numerical 
simulations in the case of strong coupling confirm the analytic 
estimates. 

PACS numbers: 72.15.Rn, 05.40.-a, 73.20.Jc 



The phenomena of Anderson localization [|] of a quan- 
tum particle in the random potential attracted a perma- 
nent interest during last 40 years. On the other hand, 
a considerable attention have been paid last years to 
the investigation of features of eigenfunctions of non- 
Hermitean Hamiltonians with disorder [p[-pT|. A popu- 
lar example of this kind is the Hatano-Nelson model , 
where in the presence of the imaginary vector-potential 
the transition between real and complex spectrum takes 
place ^,01 ■ This transition is usually referred as the 
mobility edge in Id problem (which is already surpris- 
ing iQ). In this paper we consider another example, 
where the localized(delocalized) eigenfunctions of non- 
Hermitean Hamiltonian behave in the " inverted" , com- 
pared to the Hcrmitcan case, way. 

Two recent papers deals with the non- 

Hermitean quantum mechanical Hamiltonian with an 
imaginary random potential 
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This study was motivated by the observation that the 
Euclidean evolution operator (r| exp[— ti7]|0) with H (|l|) 
after averaging over (5-correlated disordered potential V 
coincides with the probability distribution Z{r, t) for the 
Edwards self-repulsing polymer . 

Many applications of non-Hermitean operators (like 
both our examples) come from the statistical physics and 
naturally deal with the imaginary time (Euclidean) evolu- 
tion. In its turn, for the Euclidean evolution operator the 
ground state contribution (contribution from the states 
with lowest Ree) is naturally enhanced Therefore 
in this paper we are going to consider the tail of the 



eigenfunctions of the Hamiltonian (|i|), having smallest 
real part of energy Ree. In the Hermitean quantum me- 
chanics eigen-states with anomalously small e (so called 
"Lifshitz tails" |l^) originates from the rare localized 
fluctuations of the disorder V and are themselves well 
localized. For the Hamiltonian eq. (^ first of all the en- 
ergy is bounded from below Ree > The most 
surprising is the fact that for smaller and smaller Ree in 
eq. ([l]) (closer to Ree = 0) one finds (exponentially rare) 
the more and more extended states. Investigation of this 
"delocalization at the tail" in imaginary disordered po- 
tential will be the main goal of this paper jl^ . 

In all cases the tail states appear only due to some 
special rare fluctuations of the random potential. The 
standard procedure is to derive some kind of saddle-point 
equation in order to find the most probable fluctuation 
giving rise to the eigenmode with given (small) energy. 
Both the potential and the wave function in this case 
arc the instanton-like solutions of some nonlinear differ- 
ential equation |^,^. This straight scenario breaks up 
for the model eq. (1^ The necessary fluctuation in our 
case is simply V ~ const inside some sufficiently large 
connected area. The small probability to have no consid- 
erable fiuctuations on top of this V — const leads to the 
exponential suppression of density of states at the tail. 
This absence of the well defined saddle-point equation de- 
scribing the rare tail states may not be specific only for 
our Hamiltonian with random imaginary potential. It 
will be extremely interesting to find another physically 
motivated models, where the Instanton is not defined via 
the solution of any classical equation of motion. 

In this paper we will investigate in detail only the case 
of strong disorder. To this end it is natural to consider 
the discretized version of the Hamiltonian (|l|) 
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This is for 1-dimension. Generalization for the arbitrary 
dimension d is straightforward. Let random V be uni- 
formly distributed over —w/2 < Vn < w/2 with w 3> i. 
Taking into account Gaussian distributed disorder is im- 
portant if one would like to consider the mapping onto 
some kind of polymer problem. In our case however, 
Gaussian potential will only introduce unnecessary com- 
plexifications like nonuniform distribution of eigenvalues 
along the imaginary axis. In the strong coupling limit 
states in the bulk of spectrum are evidently localized, 
making it easy to distinguish the extended states at the 
border. For the continuous model (^ disorder also should 
be always considered as strong for the small enough Ree. 
Therefore although the features of solutions of the eq. (0) 
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will be very peculiar in the bulk, we still may expect some 
kind of universal behavior at the tail. 

In the leading approximation in t/w one has 
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The wave function is localized at the site m and the en- 
ergy is pure imaginary. Corrections to the energy appear 
only in the second order in t (again the generalization for 
arbitrary d is trivial) 
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There is also no real part of e. Moreover since the per- 
turbation {t) has no diagonal matrix elements the cor- 
rections to energy may be only of the even order in t 
and consequently odd in iV . Therefore, in any order of 
perturbation theory the energy is pure imaginary. The 
nonzero Ree appears if the values of potential at two 
neighboring sites occasionally turns out to be very close 
Vi — V2 ^ t. The two values of energy are now found by 
simple diagonalization of 2 x 2 matrix 
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For |Vi — V2I > 2t the energy is complex, but the per- 
turbative series in t/{Vi — V2) has finite radius of con- 
vergence. For our model (||) the Atd/w part of states are 
beyond this convergence region and acquire the nonzero 
Ree. The nonperturbative origin of Ree is clearly an 
artifact of the strong couphng limit Q. However, it 
will teach us how the extended states may appear in this 
model. On fig. 1 we have plotted the distribution of real 
part of energy for the Id and 2d versions of the Hamilto- 
nian (^) found numerically for different values w/t. 

In agreement with our prediction the huge (5-peak ap- 
pears at Ree — 0. Due to the vanishing of the derivative 
dRee/d{Vi — V2) at Vi = V2 the density of states de- 
scribed by the eq. (|^) is singular at Ree — ±t (rainbow 
singularity) 
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Here we use x and y for real and imaginary part of energy. 
The two peaks at x = ±t are clearly seen at the figure. 
(The data for smallest values of w/t corresponds more to 
the case of intermediate disorder. Still some sharp peaks 
at X close to zkt are seen in this case.) 

In the case of only two close states (^ one has — t < 
X < t. However, if where will be no disorder at all x will 
vary within —2dt < x < 2dt. These two examples show, 
how to construct the fluctuations of V which may lead 
to the smallest values of x. One should simply arrange a 
big cluster with anomalously close values of V at many 
neighboring sites of the lattice. In Id the cluster means 



a chain of q consecutive sites with close y„, « Vi+i « 
... ~ Vi+q-i ~ V. The Schrodinger equation (^) on a 
finite segment has q solutions of the form {1 < k < q) 
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Here |Kt — V| for all 1 < n < g. The tail states cor- 
respond to fc = 1. For g = 2, fc = 1 and — V2I <^t the 
eq. (^ reduces to the eq. (^). The solutions with fc > 1 
are also extended, but they have large Ree. Therefore, ipk 
with fc > 1 at given e are statistically insufficient (except 
for the highest state with k = q). Taking into account 
the mixing with sites and q + 1 may be done perturba- 
tively. However, now the correction to energy due to the 
lo(q_i_i) is pure imaginary in any order of perturbation 
theory. The only corrections to the real part of energy 
described by the eq. (Q) comes from Vn with I < n < q. 
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FIG. 1. The distribution function P{Ree) for one and two 
spatial dimensions. From top to bottom w/t = 4, 16,64 for 
d — 1 and w/t = 8, 32, 128 for d = 2. Due to the symmetry 
P{x) — P{~x) only a part of positive x region is shown for 
d = 2. J-function at x = 0, (a; ± t)~^^'^ singularity and the 
jump at a; = ±%/2i are seen at large w/t for both d = 1, 2. 

The most interesting in the eq. (|^) is that while looking 
for the state with small Ree we have got in addition the 
smooth eigenfunction extended over the large (g 3> 1) 
number of lattice sites. (We prefer not to call our states 
delocalized since the localization length ^ is usually de- 
fined via the exponent governing the decay of the tail of 
wave function and in our case ^ ^ 1.) 

On the fig. 2 we have shown the example of a rare ex- 
tended wave function for the Id model with w/t — 10. 
We have found only one such "delocalized" eigenmode 
among ~ 2 x 10^ (6 x 10^ realizations of the random po- 
tential on the sample having length / — 30 lattice sites). 
The size of the cluster is q = 11 and Ree = —1.886. 
Another quantity, which may characterize the number of 
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components of the wave function is the participation ra- 
tio p = (Y^ For the function shown 
on the fig. 2 p = 7.62, which is consistent with p — 2g/3, 
foUowing from the eq. (|^). 
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Here we keep only the leading contributions and the ar- 
gument of the exponent in the r.h.s. is found with only 
^ l/\n{w/t) accuracy. In 2d in order to get the anoma- 
lously small Ree one should have occasionally close val- 
ues of Vi within the large circle (radius r = y^q/n) . The 
analog of the eq. (||) now became 
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where jo ~ 2.4048 is the first node of the Bessel func- 
tion Jo{x). Since we are in the strong coupling regime, 
eq. (|lO|) is valid for any number of spatial dimensions. 
However, due to the another (/-dependence of e^'' ( p^ ) 
the probability now has the form (eo — —it) 



FIG. 2. The example of rare fluctuation of the potential 
V{i) (dashed) for w/t = 10, giving rise to the extended state. 
Solid line is the |V'(j)p, arbitrary units. The probability to 
find such fluctuation is P ~ 10~'. 

We devote the rest of the paper to the estimation of 
the probability to find the extended tail states. 

For any q the real part of energy described by the 
eq. (^ is bounded from below. In particular for large q 



min Ree — e, 
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In the more realistic case the real part of energy is given 
by the second order formula 



e(y)('')=4'^ + ^ 
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Here the summation goes over the states on the cluster 
(chain of the length q) described by the eq. (0) . V is the 
real potential eq. (||). The probability to find such cluster 
in the leading order may be found by simple counting 
of the number of independent variables Vn entering the 
eq. (ph. This gives 
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Here g > 1. More precisely F^?) - (x - £[,'^)(9-3)/2^ 
which gives the rainbow singularity (|6|) at a; = — t and 
a finite jump of P{x) at a; = —\/2t. For given x the 
P'''' [x) ( |lO| ) reaches maximum at q equal to 
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where eo = £\ 



{q=oo) 




'2t. The resulting probability to 



find the eigenstate, given by this largest p(*\ reads 
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In general for arbitrary d we expect (eq = —2dt) 



ln(P(a;)) In 
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On the figs. 3 and 4 we have plotted the ln(P(a;)) for 
d = 1,2, found numerically for few values of w, as a 
function of scaling variable % — [t/ {x — eo)Y^'^ . The 
agreement between the experiment and the theoretical 
formulas ( [l^ ) and (|lj) (found with only 1/ln accuracy!) 
for large w/t is rather appreciable. 

To conclude, we have considered the extended states in 
the strong coupling limit of the non-Hermitean Hamilto- 
nian with random imaginary potential. Unlike it is in the 
case of Hermitean quantum mechanics, only the extended 
states survive close to the border of spectrum (may have 
anomalously small (large) Re{e)). These extended states 
appear due to the rare fluctuations of the random poten- 
tial, where the values of V inside some sufficiently large 
area turns out to be close. To the best of our knowledge, 
this is the first consideration of the instanton-like effect, 
where the most important large fluctuation of the field 
iy) is not found by the solution of any equation of mo- 
tion. In the strong coupling limit the real part of energy 
appears only non-perturbatively. This leads to a special 
form of the density of states with various sharp features. 
This singular behavior is clearly seen at the histogram 
representing the results of numerical simulations. Sim- 
ple analytical expressions were found for the density of 
states at the tail (and correspondingly for the probability 
to find the extended states). Even though the formulas 
were found only with the logarithmic accuracy, the the- 
ory agrees well with the results of numerical experiment. 
Extension of our approach for the case of weak disorder 
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requires a further investigation. It will be also interest- 
ing to find another examples of rare "physical" events 
described by non-classical instantons. 
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FIG. 3. The logarithm of the density of states \n{P{Ree)) 
as a function of scaling variable x ~ [t/{'2t + Ree)]~^^^ for the 
Id Hamiltonian (2). From left to right w/t = 64,32, 16,8,4. 
Dashed lines are the theoretical prediction of eq. (12) with 
const = 1 and suitable overall normalization. 




FIG. 4. The same as on the fig. 3, but for 2d. Scal- 
ing variable is x = + Ree), from left to right 
w/t= 128, 64, 32, 16, 8. The theory is given by the eq. (14). 
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